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Abstract: Various implementations of the Unitarity method have been developed 
to compute one-loop amplitudes in gauge theories. In this paper we present an 
implementation which uses canonical forms to generate the rational coefficients of 
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1. Introduction 



One-loop computations are an essential ingredient in providing robust next-to-leading 
order predictions for QCD events at colliders such as the LHC [1]. A general one-loop 
amphtude for massless particles can be expressed, after an appropriate Passarino- 
Veltman reduction [2], in terms of n-point scalar integral functions with rational 
coefficients, a^, bj, Ck, 

^l-loop ^J2^J^^^J2b,^ + J2 Ck 1^2 +Rn + 0(e) . (1.1) 

The functions contain all of the logarithms and dilogarithms in the amplitude 
and Rn is the remaining rational part. The summations are over all possible integral 
functions. The Unitarity technique [3,4] determines the rational coefficients of these 
functions from the information contained in the two-particle cuts shown in fig.|l], 
using physical on-shell amphtudes as inputs. 

The cut, 

Ca,...,b = ^ j dUP^ A'''''{-ii,a,a+l, . . . ,b,i2)y<A'''''%-i2,b+l,b+2, . . . ,a-l,ii 

(1.2) 

^ where J dLIPS denotes integration over the on-shell 

phase space of the ii is equal to the leading discon- 
tinuity in the integral functions [5] of eqn. ( [L.l| ), 



b+l 



Disc 



./j a-1 \ieC' jev ke£' 

(1.3) 

TT,. 1 A i. 4.- 1 if where C, V and £' are the restricted sets of func- 
tigure 1: A two-particie cut oi ' 

a one-loop amplitude ^^^^^ ^^^^ ^ ^^^^ particular channel. In par- 

ticular, £' consists of a single term. 

The original implementation of the Unitarity method did not evaluate the cut 
directly but instead manipulated the product of tree amplitudes, using if = = 0, 
to rewrite it in the form, 



J c/LIPS^*""^(-£i, a, a + 1, . . . , 6, £2) X ^*""^(-4, 6 + 1, 6 + 2, . . . , a - 1, 



and identified the coefficients in the above with the integral coefficients of eqn. (|1.3|) . 
A feature of the original implementation is that the representation (|1.4| ) is not unique 
but one must simultaneously solve the full set of cut equations. 
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In principle, the cut momenta £i should match the momenta of the integral func- 
tions /„, i.e. they should be in 4 — 2e dimensions, however it was shown that for many 
amplitudes it is sufficient to use four dimensional tree amplitudes. Using four dimen- 
sional amplitudes allows us to evaluate the integral coefficients but not the rational 
terms Rn- For supersymmetric amplitudes, or the supersymmetric components of 
QCD amphtudes, Rn — and we term these amphtudes "cut-constructible" [3] . 

In recent years considerable advances have been made in systematising the pro- 
cess of extracting the coefficients of the basis integral functions. Progress has been 
made both via the two-particle cuts above and using generalisations of unitarity [6] 
where, for example, triple [7-10] and quadruple cuts [11] are utilised to identify the 
triangle and box coefficients. Triple and quadruple cuts are useful in that they isolate 
contributions from smaller sets of integral functions. For example, a quadruple cut 
isolates a single box coefficient. Since the cut inserts four 5-functions into the co- 
variant integral the coefficient of this box function is given by the algebraic product 
of four tree amphtudes [11] . — A one-loop amplitudes consist solely of scalar box 
functions and so quadruple cuts are sufficient to completely compute them [12-15]: 
a property shared by jV = 8 supergravity amplitudes [16]. 

The triple cut with three 5-functions is effectively a one-parameter integral which 
can be evaluated via complex methods [8,9]. One may also consider one-particle 
cuts [17]. Generalised Unitarity has been used beyond one-loop at two loops [18, 19] 
and beyond [20]. 

Our strategy is to use all possible cuts and first evaluate the box coefficients 
from quadruple cuts then the triangle coefficients from triple cuts and finally use the 
two-particle cuts to determine the bubble coefficients. This is not the only strategy 
since the two-particle cuts contain enough information to determine the coefficients 
of the box and triangle contributions as well as the bubble coefficients. The test- 
ing ground for many of these techniques has been the computation of the various 
terms in six gluon one-loop scattering amplitudes [21] both for the supersymmet- 
ric contributions [3,4,22,23] and also for the "cut-constructible" parts of the QCD 
amplitudes [4,24,25]. For full QCD amplitudes the rational terms must also be cal- 
culated. Unitarity can be used to determine these, however this requires the use of 
tree amplitudes defined in D = 4 — 2e dimensions [26] . Alternatively, the rational 
pieces can be obtained using on-shell recursion [27], a method akin to that for tree 
amplitudes [28]. This has a numerical implementation together with Generalised 
Unitarity [29] . Alternate numerical implementations exist for variants of this strat- 
egy [30]. One may also use specialised Feynman diagram techniques which focus on 
the rational terms [31-33] . 

The approach we adopt recognises that there are a limited number of distinct 
structures that appear in the cut integrals. These may be evaluated in a number of 
ways: conventional covariant integration, fermionic integration [23,25], direct extrac- 
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tion [8] or integrand level reduction [34]. By determining the contribution of each 
structure to the relevant coefficients we construct a canonical basis which, once con- 
structed, can be used for any amplitude. While this approach reproduces results from 
other methods, the decomposition into canonical forms is carried out directly on the 
four dimensional tree amplitudes without re-parametrising the dLIPS integration. 
It produces compact, explicitly rational expressions for the integral coefficients. We 
illustrate this process by presenting the J\f — 1 contribution to one-loop seven gluon 
scattering in closed analytic form^. 



2. Organisation of the Amplitudes 

The organisation of loop amplitudes into physical sub-amplitudes is an important 
step toward computing these amplitudes: although eventually all the pieces must 
be reassembled. For one-loop amplitudes with adjoint particles, one may perform 
a colour decomposition similar to the tree- level decomposition [35]. This one-loop 
decomposition is [36], 

Ln./2J+1 

^ J2 Gr„;e (<7) ^„^,(<7) , (2.1) 

C=l aeSn/Sn;c 

where [x\ is the largest integer less than or equal to x.^ The leading colour-structure 
factor, 

Gr„;i(l) = N, Tr(r'^i---r'^") , (2.2) 

is just Nc times the tree colour factor and the sub-leading colour structures (c > 1) 
are given by, 

Gr„.e(l) = Tr (T"i • • • T"=-i) Tr (r"= • • • T"") . (2.3) 

Sn is the set of all permutations of n objects and Sn;c is the subset leaving Gr„.c 
invariant [36] . The contributions from fundamental representation quarks circulating 
in the loop can be obtained from the same partial amplitudes, except that the sum 
only runs over the An-i and the overall factor of in Gr„.i is dropped. For one-loop 
amplitudes of gluons the An-c, c > 1 can be obtained from the An-^i by summing over 
permutations [36,3]. Hence it is sufficient to compute An-i in what follows and, for 
clarity, we refer to these as The partial amplitudes /!„ have cyclic symmetry 
rather than full crossing symmetry. 

We choose to use a super symmetric decomposition. Instead of calculating the 
one-loop contributions from massless gluons, A^n , or quarks, A^n"^^ , circulating in the 

^These are available in Mathematica format at http://pyweb.swan.ac. uk/'^dunbar/sevengluon.html 
^We have inserted a factor of i in this definition to avoid universal factors of i appearing in our 
explicit formulae. 
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loop, it is considerably more convenient to calculate the contributions from the full 
A/" = 4 multiplet, a A/" = 1 chiral multiplet and a complex scalar circulating in the 
loop. In terms of these, 

/I [1/2] _ J Ar=l chiral _ M ^ ' > 

— ■ 

The amphtudes are also organised according to the helicities of the outgoing 
gluons which may be ±. We use polarisation tensors formed from Weyl spinors [37], 

where k is the gluon momentum and q is an arbitrary null 'reference momentum' 
which drops out of the final gauge-invariant amplitudes. The plus and minus labels 
on the polarization vectors refer to the gluon helicities and we use the notation (ij) = 
{k~\k'j) , [ij] = {kf\kj). In twistor-inspired studies of gauge theory amphtudes [38] 
the two component Weyl spinors are often expressed as, 

A„=|^+), \a=\k-). (2.6) 

Helicity amphtudes are related to those with all legs of opposite helicity by conjuga- 
tion: [ah) [6a]. Using spinor helicity leads to amplitudes which arc functions of 
the spinor variables {ah) and [ah] and combinations such as {k'l\'^\k'^) — [ki\p\kj) = 
[ki'p\{pkj). It is useful to define combinations of spinor products, 

[a\P^...f\m) = [ab]{bm) + ■■■ + [a f]{f m) . (2.7) 

In this article we complete the computation of the A/" = 1 contributions to 
seven gluon scattering. Up to conjugation and relabeling, there are nine inde- 
pendent helicity configurations for the colour ordered amplitudes. The amplitudes 
^7(1+, 2+, 3+, 4+, 5+, 6+, 7+) and ^7(1", 2+, 3+, 4+, 5+, 6+, 7+) vanish to all orders in 
perturbation theory within any supersymmetry theory so, 

A^=\l+, 2+, 3+, 4+, 5+, 6+, 7+) = A^=\l- , 2+, 3+, 4+, 5+, 6+, 7+) = . (2.8) 

Amplitudes with exactly two negative helicities are referred to as MHV ( "maximally 

helicity violating") amplitudes and those with three negative helicities as NMHV 
( "next to MHV" ) amplitudes. There are three independent MHV and four indepen- 
dent NMHV helicity configurations for the seven gluon amplitude. The seven-point 
MHV amplitudes [4] and the NMHV amplitude with the three negative helicity legs 
adjacent =^(1", 2", 3", 4+, 5+, 6+, 7+) [7] are specific cases of known all-n expres- 
sions. We present explicit forms for the remaining three NMHV partial amplitudes 
which will be made available at http://pyweb.swan.ac.uk/'^dunbar/sevengluon.html 
in Mathematica format. 
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3. Canonical Basis for Bubble Coefficients from Two-Particle 
Cuts 



d^£5if)5ii£ - Pf) ^^ y (3.2) 



Consider the general decomposition of the product of tree amphtudes appearing in 
a two-particle cut in terms of canonical forms 

^^■■^^(-^i,---,^) X A*'-^'=(-4,---,^i) = Y^c^m^)^ (3-1) 

where the Cj are coefficients independent of The forms jFj must have zero spinor 
weight in the i.e. they must be invariant under \£j) — e^'^^\£j), \£j\ — > e~^'^^\£j\. 
Each two-particle cut receives contributions from box functions, triangle functions 
and a bubble. The coefficients of the box functions are most simply determined from 
quadruple cuts and the triangle coefficients from triple cuts, so we only use the two- 
particle cut to determine the bubble coefficients. Consequently, we must determine 
the contributions to the bubble coefficients from the various canonical forms, !Fi. 

We examine the !Fi according to their leading order in the cut momenta £. For 
a generic QCD amplitude the contributions have a maximum order of In the 

A/" = 1 contributions to Yang-Mills amplitudes cancellations generically reduce this 
to order . Individual terms of order with < give no contribution to the 
bubble coefficients. To see this, we manipulate the cut into terms of the form, 

P^+«(£) 

where the Qi may be null or non-null. For iV < 0, Passcrino-Veltman reduction on the 
corresponding covariant integral yields box and triangle integrals only. Consequently 
we are only interested in terms of order with N > 0. 

3.1 Order ^° Terms 

These are the canonical forms which are required to obtain the J\f = 1 contributions 
to Yang-Mills amplitudes. 

Functionally we start with the simplest non-trivial case, 

where kA is taken to be real. There are many ways to evaluate the contribution to 
the bubble coefficient from this simple form. We chose to manipulate this as if it 
were a covariant integral. This means effectively replacing, 

dLIPS ^ y ^ , (3-4) 

then evaluating the covariant integral only keeping the coefficient of — ln(— P-^) in 
the result. The integral is then a hnear triangle integral with massless leg kA as 
shown. 
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kA 

which can be evaluated by shifting the momenta using Feynman parameters, 



'^^i';- fc>2 + P'^as, (3.5) 



where P = ka + ka+i + ■ ■ ■ k^ is the total momentum across the cut and as refers to 
the Feynman parameter of the £2 propagator. 

The a2 term drops out of the integral and we use, 

. H^n-H-iP~>^^X ,3,, 

The coefficient of — ln(— P^) is the bubble coefficient and thus Tii evaluates to Hi, 

Note we assume in the above that /ca is real. If A denoted a complex combination of 
momenta then 1^4] should be replaced by 1^4)* in the canonical form. We also have 
the conjugate result, 

H.(C;A«^|§f|^/^dC;fl;P|^|^. (3.8) 



These forms satisfy. 



Hi[ P\A]; B- P] = -^HM- P\B)-P]. (3.9) 



We have chosen to determine the contribution from this form using covariant 
integrals, of course one obtains the same result by applying fermionic integration [23] 
or direct evaluation [8]. Our aim is to use this simple result for Hi as the starting 
point for generating many further terms. First consider 7i to be a holomorphic 
function of one of the £i, i.e. H, = TC{\i)), then we can define, 

By spinor weight the number of Ai is the number of Bj. The cases where mul- 
tiple poles appear are treated separately: such terms cancel in amplitudes by the 
factorisation theorems. 



By expressing this as a partial fraction we can split Hn algebraically into a sum 
of Til terms, 

nn{A;B,;i) = = ^c,?ii(A;i^i;^), (3.11) 



where the coefficients q are given by, 

_ n?=2(^.A) 

The bubble coefficient generated by Hn is thus. 



(3.12) 



HnlA, B^, P] - ^ n,,.(^.A) {A\P\A.Y ^'-''^ 

The same formula applies whether we have a holomorphic expression in ii or £2- 
When we have a mixed expression. 



nr=i(-B.^i) uT=i{c,i2) uuiEkii] uu[Gii2] 
nr=i(A^i) n7=i(^.4) nLJ^^^i niAHihy 

we first rewrite it as a holomorphic expression in £1 and £2, 

nr=i(A^i) n;ii(^.^2) rlj^^i^k^) nij^d^Ki) 

then use the identity, 

(ah) (ah) P'^iab) 



(3.14) 



{bi2) {bii) {bii)[ii\P\b) 

to replace it by. 



(3.15) 



(3.16) 



nLi(A4) U%AD,^.) ULAF^W^) UlAHi\P\i.) ^'''"^ • ^'-''^ 

Only the leading term contributes to the bubble coefficient and it gives an overall 
contribution of Hn+m+p+q[A: Dj, P\Fk\, P\Hi]; Bi,Cj, P\Ek\, P\Gi\; P], provided the 
Ai, Dj, and P\Hi] are all distinct. 

We can have terms with D — A. In this case we decompose into terms of the 
form. 

This special case canonical form gives, 

H|KAi^.,i^.P]^ '^l™7-> . (3.19) 
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There is a second class of functions arising from terms with propagators involving 
non-null momenta such as, 

go{B-D-Q-h) = (3.20) 
where ^ q. We can relate this to the Tii form using the identity, 



leading to, 

1 |B|,.|B) = _S^Ml£M+,„b-leadmg. (3.22) 



We then make the replacement, 

{i\PQ\i) ~ {^\PQ\(.) = {eP)[PQ]{Qi), (3.23) 
where P and Q are the null linear combinations of P and Q, 

= -^(pV-(i^-Q-^)P^), Q'' = ^j={-P'Q''+{P-Q+^)P^), 



(3.24) 

with A3 = AsiP, Q) = 4(P ■ Qf - AP'^Q^. A3 is the Gram determinant of the three- 
mass triangle integral having legs of momenta P, Q and — P — Q. The leading term 
in eqn. (|3.22| ) then has precisely the form of an 7^2 function. Splitting this into a pair 
of Hi functions gives two terms that are not individually rational because P and Q 
contain factors of y/A^. The pair of terms are however the irrational conjugates of 
each other, so the sum is rational. To have canonical forms which yield explicitly 
rational coefficients we choose to evaluate as a separate canonical form which is 
manifestly rational: 

A3 A3 

We commonly find the form, 

GM. Bo, Pi, P, Q, ^1) - ^^^^ . (3.26) 

For {AP), {AQ) ^ 0, this can be decomposed as an function but at the cost of 
introducing irrational coefficients. Once again, combining these terms generates a 
manifestly rational form. For {AP), (AQ) ^ 0, 

[P|P[P,Q]|A)(Pi|[P,Q]|Po) 



Gi[A;Po,5i;P;Q;P] = 

+ 



2A3(A|pg|A) 

[D\P\A){{BoA){Bi\P\A] + {BiA){Bo\P\A]) ^^'^^^ 
2{A\PQ\A){A\P\A] 
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We can extend this form to determine the bubble contributions arising from, 

G(A-B H-n-n-f)- ^ mi\Bo)ir:=i{^^B,) 

g^iA, Bo, B,, D, Q, 4) - ^^-^ — nr=i(^iA) — ' ^ ^ ^ 

by sphtting it into a sum of Qi terms, just as we spht Hn into a sum of Hi terms, 
GniA; Bo, Bi, D; Q; £i) = ^ Q GiiA; Bo, S„; D; Q; (3.29) 



with. 



The Hn and Gn functions are sufficient to evaluate the bubble coefficients of the 
Af = 1 contributions to Yang-Mills amplitudes. We will show this by example in 
the following sections where we explicitly evaluate the seven-point J\f = 1 one-loop 
contributions. 

In general we may also have terms with multiple propagators involving non-null 
momenta, 

m 



Using the constraint (^ — P)^ = 0, we have. 



(3.31) 



{e + Qif = [i\Qi + ^P\i) ^ mi\i), (3.32) 

where Qi is a non-null linear combination of Qi and P. This allows any multiple 
propagator terms to be written as, 

mi\mQ2\i)..- ^ ■ ^ 

Partial fractioning on the A(£)'s then gives a sum of terms of the form, 

m 

mi\i)miQ2\i)...^ [^|Qi|^)(^Qi)[Qi4](Q2^)...' 



(3.34) 



which can be further split into a sum of Qi forms at the expense of introducing 
irrational factors in ^/A^(Q~'Q~). As with the ■\/A3(P, Q) terms, these all arise in 
irrational conjugate pairs and the sum is rational. 

3.2 Terms of Order and i'^ 

For the scalar contributions to Yang-Mills amphtudes we need forms of order and 
i'^. In general we will needs forms denoted for contributions of order where 
the denominator has n factors of ( £) . 



- 10 - 



The higher order Ti, and H forms are: 



nl = [D\£,\B) ^ H'o[B;D;P] = l[D\P\B), 
^,_ [D\ii\B,){e,B2) , 

El, B,- D- P] = {[D\A\B,)[A\P\B,) + [D\A\B,)[A\P\B,)) 

+ 4[A|P|A) (iD\P\B^)[A\P\B,) + [D\P\B,)[A\P\B,)) , 

Hl[Br,B^;D^,D2,P] = ^[A|P|Si)[L>2|P|S2) + ^[A L'2](Si 

Hi = 



{iiA) 



Hl[A; Bi, B„ B,- D„ D,; P] = ^^^^^^^^^^3 i[D,\A\B,)[D2\A\B,)[A\P\B,) + V,iB,)) 
+ oJ!'2.^2 m\P\Bi)mA\B,)[A\P\B,) + Pi2(5^, A)) 



2>&[A\P\AY 
1 

IS[A\P\A) 



+ .or.lOlM ([A|P|Si)[L>2|P|S2)[A|P|S3) +P6(S0) , 



(3.35) 

where Vq{B^ and Vi2{Bi,Di) represent a total of six and twelve permutations re- 
spectively. 

The higher order terms involving propagators with non-null momenta can be 
reduced to or Ql forms, where, 

^ ^ m\i){Ae) ' ^'^•'^^^ 

and is a polynomial of degree n in i. To evaluate these forms we again make 

use of the identity: 

[Di] [D\Q\^) [D\m _ (3_37) 



m\i) [e\Q\i){i\PQ\i) {i\PQ\i)' 
This allows us to write, 
g^{A\PQ\A) 

_ P'fi£){B£) ( [D\Q\A){CA) 4 / [D\Q\P){C\PQ\A) [D\Q\Q) {C\QP\A) \ 

,.,wp,. 4 ( [D\P\P){C\PQ\A) [D\P\Q){C\QP\A) \ 

(3.38) 
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This splits the term of interest into pieces with the same overall power count in d but 
only masslcss propagators and terms with a reduced power count in I involving the 
original propagator. The former are readily evaluated using our basis of H functions, 
while the latter rely on G functions. The terms involving the propagator along with 
{P£)~^ or {Q£)~^ give rise to special cases of the G functions as discussed below. 
Using this splitting procedure, the order £^ term, Ql, is expressed in terms of a sum 
involving special cases of the order £° Qi form, namely, and Q[, where, 

^ [Di]{CmBi) 

m\e)m 

Evaluating this form gives the bubble contribution, 

_^ mQ,mC)[Q\P\B) + {P'[D\Q\B) - {RQ + ,/A;/2)[D\P\B))[Q\P\C) 

P^s • 

(3.40) 

Q[ and Gf are obtained by irrational conjugation. 

Similarly for the order terms we require special cases of Ql with {D Q) = 
and (DP) = 0. Setting f(£) = [E\£\F), we have, 

,4 ^ [E\£\F)[D\£\C){£B) 

[£\Q\£){£Q) 

Making the definitions, 

/3 = [E\P\F)[D\P\C)[Q\P\B) - [E\Q\F)[D\Q\C)[Q\P\B) 

-[E\Q\B)[D\Q\C)[Q\P\F) - [D\Q\B)[E\Q\F)[Q\P\C) 

and 



(3.42) 



(3.43) 



7 = [{[E\P\F) - [E\Q\F)) {[D\P\C) - [D\Q\C))[Q\P\B) 

- {[E\P\C)[D\Q\F) + [E\Q\C)[D\P\F))[Q\P\B) 
-[E\Q\B)[[D\P\C)-[D\Q\C))[Q\P\F) 

-[D\Q\B)([E\P\F)-[E\Q\F)yQ\P\C) 
we find that this form gives a contribution to the bubble coefficient of, 

Gr« = -8^5^±pi. (3.44) 

Again Ql'^ and G];^ are obtained by irrational conjugation. 

Although these forms contain -\/A^, we are always interested in the sum of 
irrational-conjugate pairs and the final canonical form is guaranteed to be rational. 
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4. Example: Seven-point M =1 Contributions 



The basic NMHV amplitudes are: 

A:A(r,2-,3+,4-,5+,6+,7+) 
i?:A(l-,2-,3+,4+,5-,6+,7+) 
C:A7(l-,2+,3-,4+,5-,6+,7+) 
Z}:A7(l-,2~,3-,4+,5+,6+,7+) 

Amplitude A has no permutation symmetries while amplitudes B and C have the 
following invariances: 

B : (1234567) ^ (2176543) 

C : (1234567) ^ (5432176) 

Amplitude D has the simplest structure and is an example of a "split-helicity" am- 
plitude whose n-point expression is given in [24]. 

The supersymmetric cancellations present in a A/" = 1 computation lead to the 
one-loop amplitude being of the form, 

= E «^ ^4 + E ^3 + E ^2 , (4.2) 

with no further rational terms. Since all the coefficients can be evaluated from four- 
dimensional unitarity these contributions are termed "cut constructible" . 

For the contribution from the A/" = 1 chiral multiplet the coefficients of the box 
integral functions contain sufficient information to determine the coefficients of the 
one and two mass triangle functions. As discussed in the Appendix we choose to 
absorb these triangles into the box functions, leaving a basis of truncated boxes J^l, 
three mass triangles and bubble functions, 

ieC ieDg™ ke£ 

4.1 A^=i(l-,2-,3+,4-,5+,6+,7+) 

This amplitude can be decomposed into 20 integral functions: 

/iA/'=l chiral/1- 9- 0+ 4- cr+ ^+ ^+\ _ A -pSm 
^7 l-L , ^ 5 ^ , , D 5 ' J — "1 •^6{71}{23}{45} 

I A 'r'2m h , A r2m h , A r2m h , A ^2m h 

"T "-2 1{23}{456}7 -T a^J- 3{45}{671}2 "T "4 ^ 3{456}{71}2 T" "5 5{671}{23}4 

I ^^-r-2me _|_^-4-r-lm _|_ -r-lm 

^ "6 3{45}6{712} ^ "7 '^234{5671} ^ "S 345{6712} 

I lA rSm tri , lA jSm tri 

+ '^1 -'{23}{45}{671} + '^2 -'{71}{23}{456} 

+ cfhitus) + 4h{t2M) + 0^/2(^345) + cf /2(t456) + 0^/2(^671) 
+ C^/2(t712) 

+ 4h{s2z) + 4^2(534) + df 12(345) + dfhisn), 

(4.4) 
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where we have chosen to label the boxes and three-mass triangles by the clustering 
of the legs. 



4.2 A^=i(l~,2-,3+,4+,5-,6+,7+) 

This amplitude can be decomposed into 25 integral functions: 

^7 U , ^ ,-3 , ^ , ^) , , ' J — ^1 -^6{71}{23}{45} + ^2 '^4{56}{71}{23} 

I B -r2m h I B -r2m h < B •r2m h < B •r2m h 

"T "3 1{23}{456}7 "T "4 3{45}{671}2 "5 6{71}{234}5 "6 1{234}{56}7 

I B ■r2m h I B 'r2m h 

"T '^T 5{671}{23}4 "s 3{456}{71}2 

I B ■r2m e , B r2m e , B -rim 

"T Ut^J' 3{45}6{712} "T "10-^ 4{56}7{123} "ll-^ 456{7123} 

I rSm tri 1 t_B T3m tri 1 l_B 7'3m tri 

+ "1 -'{23}{45}{671} + "2 -'{71}{23}{456} + "3 -'{56}{71}{234} 

+ cfhitus) + 0^/2(^234) + cf 12(^345) + cf /2(t456) + C^hihej) 
+ C^f /2(S23) + C?f /2(S45) + Hs^q) + C^f /2(s7l). 

(4.5) 

4.3 A^=i(l-,2+,3-,4+,5-,6+,7+) 

This amplitude can be decomposed into 37 integral functions: 

^7 U , Z , J , 4 , , , ^ J — «1 •>^6{71}{23}{45} + "2 •^7{12}{34}{56} 

I C -r2m /i I C -r2m /i 1 C -r2rra /i _i_ C r2rn h _i_ ■T'2m h 

-h ag ^i{23}{456}7 "T "4 2{34}{567}1 3{45}{671}2 "6 •^4{56}{712}3 "7 6{71}{234}5 

I C -772m /i I C -772m /i I C 'r2m h , C r2m h , C -772m h 

"T "8 1{234}{56}7 "T "9 3{456}{71}2 "lO-' 4{567}{12}3 "ll-^ 5{671}{23}4 "12-' 6{712}{34}5 

I C ^2m e _L /i'-^ -772m e 

"T "l3-^ 4{56}7{123} "M-' 6{71}2{345} 

I C ^Im I C ^Im I C ^Im , C ^Im , C ^Im 

"T "is-' 123{4567} "T "le-' 234{5671} "17-' 345{6712} "r "iS-' 456{7123} "iQ-' 712{3456} 

C rSm tri , r C 7-3™, tri , iC rSm tri , iC rSm tri , iC jSm tri 



I to T6m tri I r O jdm tri 1 7 O T6m tri 1 r O rdm tri 1 r O rJm trj 

+ "1 -'{12}{34}{567} + "2 -'{712}{34}{56} + '^3 {71}{23}{456} + -'{71}{234}{56} + "5 -'{671}{23}{45} 

^567) 



+ Cf /2(tl23) + C^h{t2M) + 0^/2(^345) + C^ihiUm) + C^hih 



+ C^/2(t67l) +0^/2(^712) 

+ d^h{si2) + d^h{s23) + rf^/2(s34) + rff /2(s45) + d^ his^o) + d^hisn). 

(4.6) 

Altogether these NMHV A/" = 1 contributions are specified by 82 coefficients. 
The box and three-mass triangle coefficients are special cases of generic forms which 
are given in appendix |C| and section ^ respectively. The twenty Cj coefficients are 
either special cases of the generic Cq function ( |C.7|) or are given by one of four forms, 
Cx, specific to the seven-point case. The remaining parts of the amplitudes are given 
by the fourteen df functions, which are not all independent but can be expressed in 
terms of six Dx functions. 
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4.4 Cx Functions 

We illustrate the calculation of the Cx functions by considering an explicit realisation 
of the Cb function which arises as the coefficient of 12(^712 ) and is obtained by 
computing the ^712 cut of ^47(1, 2, 3~, 4+, 5~, 6+, 7), where precisely one of legs 7, 1 
or 2 has negative helicity. We label this leg mi. 

The product of tree amplitudes generated by the cut is, 

^A*-<'(-£^7,l,2,£2') X A*-<=(-4,3-,4+,5-,6+,£r'), (4.7) 
h 

where the summation is over the complex scalar and fermionic states of the M = 1 
chiral multiplet. The six point NMHV tree amplitude has three terms, 

3^, 4+, 5-, 6+, ^r'') = + T!^ + , (4.8) 

where. 





|P456|5)'[^l|P456|5)' 


t456[^i4][43](4 5)(5 6)[£i| 


|A56|4)[3|P456|6) 



[^2|P456|5)''' 



[6|P56£j3)^((£2 3)[6£i])^ / [6|P56.j3) y 

' t56£,(^23)(34)[56][64][5|P56£j^2)[4|P56^i|4) V {i2S)[Q£i]J ' 



T.: 



[4|P345|^l)'[4|P345|^2)=^ / " [4| P345 



t345(6£i)(£i4)[34][45][3|P345|6)[5|P345|4) V [4IP345K2) 



h 



with h denoting the helicity of the leg £2: h = for a scalar and /i = ±1 for a 
fermion. 

Summing over the multiplet in eqn. ( |4.7| ) leads to term by term cancellations 
which we can express as p factors multiplying each product of h = terms, leading 
to a cut integrand of the form. 



(mi£i)2(mi4) 



(71) (12) (2 £2) (^2^ 
where the p factors are. 



;^x(T0p,+T2V2 + T0p3). 



(4.10) 



Pi 



(mi|P7i2P456|5)^ 



(mi£2)[4|P456|5)(mi£i)[£l|P456|5)' 



(mi£2)[6|P56£j3)(mi£i)(£23)[6£i]' ^ ' ^ 



P3 

with. 



(mi£i)[4|P345|^l)(mi£2)[4|P345|4)' 



|>B2) = [6 5](5 3)|mi) + [6 7](mi3)|7) + [61](mi3)|l) + [6 2] (mi 3)|2). (4.12) 
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Consequently the contribution of the Ti terms to the cut integrand is, 
{mi£i){mi £2) 



X 

\2 





|P456|5)[^l| 


|P456|5)(mi|Pn2P456|5)' 


i456[^i4][43](4 5)(5 6)[^i| 


|i'456|4)[3|P456|6) 



(mi|Pn2P456|5)^ ^ . .... . _, ^ 

X ; ~ — ; — ~ ; X 



(7 2) (4 5) (5 6) [3 1 P4561 6)^456^712 (2 £2) {ii 7) 

(mi|Pn2P456|5)2 ^ (mi £1) (mi £2) ^ 



[41 


\P456\5)[il\ 


-P456 5) 


[£2 3] [41 


1^45614) 



(41 


-P712-P456 5)(£2 


-P712-P456 5) 


(41 


|P712|3](£2| 


IP712P456I4) 



(7 1) (1 2) (4 5) (5 6) [3|P456|6)t456t7i2 (2 £2) {£1 7) 

(4.13) 

This can be recognised as a Ti^ canonical form and thus yields a contribution to the 
bubble coefficient of, 

(mi|Pn2P456|5)=^ 



(71)(12)(45)(56)[3|P456|6)t 4564l2 

X H^[2, 7, Pvi2|3], P712P456I4); mi, mi, P712P456I5), P712P456I5); P712]. 

(4.14) 

Similarly the T2 and T3 terms give a contribution, 

6-4(2, 7, P34I5], P7i2n6|4); mi, mi, 3, Yb2, Yb2] 6; P34; P712] 
H^[2, 7, 6, P345I5]; mi, mi, P345I4], P345I4]; P7i2]. 



(71) (12) (3 4) [5 6] 

[4|P456|mi)2 



(71)(12)[34][45][3|P345|6)i345' 

(4.15) 

We define the Cb function by generalising the result of this specific cut: 
CB[a,6,c,d,e,/,5r;mi] = 



(mi 


PgabPdef \ 


\e? 


{ga){ah){de){ef)[c\ 


\P 


'def\f)tdeftgab 



{ga){ab){cd)[ef] 



X H^lb, g, Pgab\c], PgabPdef \d); mi, mi, PgabPdef\e), PgabPdef Pgab] 

Gilb, g, Pcd\e\, PgabPef\d);mi, mi, c, Yb2, >s2; /; Pcd, Pgab] 



[d\Pcde\ 


mi)^ 


{g a){ab)[cd][de][c\ 


\Pcde\f)^cde 



Hi[b,g,f, Pcde I e] ; mi , mi , Pcde \ d] , Pcde \ d] ; Pgab] , 

(4.16) 
with, 

\Yb2) = [/e](ec)|mi) + [/^](mic)|^) + [fa]{rmc)\a) + [f b]{m,c)\b). (4.17) 

Five of the bubble coefficients can be expressed in terms of the Cb function: 

= Cb[1, 2, 3, 4, 5, 6, 7; 1], = -Cb[7, 6, 5, 4, 3, 2, 1; 1], 
= Cb[6, 7, 1, 2, 3, 4, 5; 5], = -Cb[5, 4, 3, 2, 1, 7, 6; 5], (4.18) 
c^ = Cb[7,1,2,3,4,5,6; 1]. 
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We define three further functions in this class, Ca, Cc and Cd] 



CA[a,b,c,d,e,f,g;mi] = 

{mi\PgabPcd\ey 



(c d){de)[f\Pd„\c){ga){a h) tcdetgab 
[d\Pef\m,)' 



H3 [b, g , Pyab I /] ; mi , mi , PgabPcd I e) ; Pgab] 



(c/ a) (a 6) e] [e /] [/ 1 Pde I c)tde/ 

(4.19) 



Cc[a,6, c, rf, e, f,g;mi\ = 
[def{rm ff 



tcde[cd] [c\Pde\f){ga){ah) 

{m^\PgabPde\f? 



Hsig, b, Pcd\e\; /, mi, mi; Pgab] 



+ 



tde ft gab {de){ef)[c\Pde\f){ga){ab) 

X Hi[c, Pef\d), Pga\b), Pablg); Pde\f): Pde\ f) : Pgab\mi) , Pgab\mi) ] Pgab] 

1 



[e f]{cd) {g a) {a b) 

X G^[Pcd\e],PgabPef\d)^ -Pgafel/], 9^ c, mi, mi, Pgab\e],Yc, Yc\ e; Pcd\ Pgab], 

(4.20) 



CD[a,b,c,d,e, f,g;mi] = 

{dey[f\Pgab\m^)^ 

tcdetgab{cd) [f\Pde\c) {g a) {tt b) 

[f\Pde\m,y 



H3[b,g, PgabPcd I e) ; mi , mi , Pgab \ f] ; Pgab] 
H4[b,c,g, Pef I d] ; mi , mi , P^e I /] , Pde I /] ; Pgab] 



tdef[d e] [e f][f\Pde\c){ga){ab)' 

-G5[f, b, g, Pef\d], PgabPcd\e); e, e, mi, mi. Yd, Yd; c; Pef, Pgab], 



{ef)[cd]{ga){ab) 

(4.21) 

where. 



\Yc) =[e/](/c)|mi) + (mi c){[eg]\g) + [ea]\a) + [eb]\b)), 
\Yd) = - [cd]{de)\m,) + {em,)([cg]\g) + [ca]\a) + [cb]\b)). 
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The remaining cf coefficients are then: 

cf = Co[2, 3, 4, 5, 6, 7, 1; 3, 4], = Co[3, 4, 5, 6, 7, 1, 2; 3, 1], 
4 = Ca[4, 5, 6, 7, 1, 2, 3; 4], = C^[5, 6, 7, 1, 2, 3, 4; 4], 

c^=Co[l,2,3,4,5,6,7; 7,4], 

cf = Co[2, 3, 4, 5, 6, 7, 1; 3, 5], cf = Cc[3, 4, 5, 6, 7, 1, 2; 2], 

cf = C^[4,5,6,7,1,2,3;5], cf = C,3[5, 6, 7, 1, 2, 3, 4; 5], (4.23) 
ci = -Ca[Q, 5, 4, 3, 2, 1, 7; 5], cf = -Cc[7, 6, 5, 4, 3, 2, 1; 1], 

cf = -Co[l,7,6,5,4,3,2;7,5], 

cf = Co[2, 3, 4, 5, 6, 7, 1; 5, 2], cf = Cc[3, 4, 5, 6, 7, 1, 2; 3], 

C3^=Co[4,5,6,7,l,2,3;4, 1]. 

4.5 -Dx Functions 

The Dx functions arise as the coefficients of l2{sab) when we consider cuts of the 
form, 

J dUPS A''^{-iu a, b, i2) X A''^{-i2, c, d, e, /, g, ii) D[a, b, c, d, e, /, g] . (4.24) 

In order to obtain non- vanishing A/" = 1 and scalar contributions, legs a and b must 
be of opposite helicity. For the J\f = 1 contribution the two helicity configurations 
for a and b are trivially related: 

D[a'^,b'',c,d,eJ,g] = -D[b' , , c, d, ej, g] . (4.25) 

Consequently, the number of independent Dx functions corresponds to the number of 
independent helicity configurations for the legs c, d, e, /, g that contain two negative 

and three positive helicities. There are six such configurations: ( h + +), ( — h 

- + +),(- + + -+),(- + + + -), (+-- + +) and (+ - + - +). 

To evaluate these cuts we need the explicit forms of the seven-point NMHV tree 
amplitudes where two external states are scalars or fermions (given explicitly in Ap- 
pendix |D|). Using these forms for the tree amplitude A^'^^^{i2, c~, d~, e"*", /"*", g~^, ii^), 
we express the cut of the ffist of these helicity configurations in terms of canonical 
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forms and obtain, 
DA[a,h, c,d,ej, g] 

-H2[b,g;a,Pcde\e\;Pab\ 

ifd) 



[e\ 


Pcde 


|a)2 


{ab){fg)[cd][de][c\ 


Pcde 1 f ) ^cde 



_ [ah]{de){ef){fg)[c\P,,\f) 

Hi[a, C, Pcde I/), PabPabcPde\f);b, g, Xa, Xa] Pab] 

1 (4.26) 

2_ [ah]{de){ef){fg)[c\Pae\f) 

H, [a, Pab\g) , C, Pcdel /) , PabPabcPde \ f); b, Pab\ f) , Pef M) , X^, X^; Pab] 

{cdYighf 



[ah]{cd) {de){e f)[g\Pab\c)tgab 
{ac^glPaM^ 



{ab){de){ef)[g\Pab\c)[g\ Pabc I d) tabctdef 
where, 



H2[b,PabPgab\f);a,Pab\g];Pab] 

H2 [b, PabcPde \f)]a,C] Pab] 



\Xa] ={fd)[ba]{ag)\g] + {f d)[ba]{a f)\f] - {a f)[ba]{de)\e] 

+ i{f9)[9c]{cd) + [c\Pabc\f){cd) + {fd)Sab)\b]. 



(4.27) 



For the other hehcity configurations we define Db,c,d,e,f in a similar fashion. 
The exphcit forms of these are given in appendix ^ In terms of these we have, 

4 = Dd[2, 3, 4, 5, 6, 7, 1], df = -Da[A, 3, 2, 1, 7, 6, 5], 
df = -De[5, 4, 3, 2, 1, 7, 6], c?^ = Db[7, 1, 2, 3, 4, 5, 6], 

d^ = -Dc[2, 1, 7, 6, 5, 4, 3], d^ = De[A, 5, 3, 2, 1, 7, 6], 

rff = -De[Q, 5, 7, 1, 2, 3, 4], = Dc7[l, 2, 3, 4, 5, 6, 7], (4.28) 

df = -Db[2, 1, 3, 4, 5, 6, 7], d^ = Dc[3, 1, 7, 6, 5, 4, 2], 
= -Dc[3, 5, 6, 7, 1, 2, 4], = ^^[4, 5, 3, 2, 1, 7, 6], 
d^ = Df[5, 6, 7, 1, 2, 3, 4], d^ = -Df[1, 7, 6, 5, 4, 3, 2]. 

For each amphtude we have checked at exphcit kinematic points that these bub- 
ble coefficients satisfy, 

^ Ci + ^ rfj = (4.29) 

i 3 

These conditions ensure that each amplitude has the correct 1/e IR singularity [39]. 
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5. Canonical Basis for Triangle Coefficients from Triple Cuts 



Generalisations of unitarity can be used to determine the coefficients of triangle and 
box functions simply. If we consider a triple cut, 

C3 = - / dHS{il)S{il)S{il)Y,A(i-io),mr--,j-l,{ii)) 

J ^ ' (5.1) 

X a((-£i),j, 1,(^2)) X ^((-£2), m-l,(£o)), 

The minus sign in this equation is for when using colour-ordered partial amplitudes 



as normalised in eqn. ( |2.1| ). This has contributions from the discontinuities of a single 



triangle and several box functions, 




Disc 



(5.2) 



The information in this cut may be used to determine the triangle coefficient hy 

One of the advantages of the supersymmetric decomposition of gluonic ampli- 
tudes is that the one and two-mass triangle coefficients need not be explicitly com- 
puted: for the A/" = 4 contributions they are absent while for the M = \ and scalar 
contributions they are tied to the box coefficients and can be incorporated into the 
truncated box functions. Consequently this section will focus on the case where 
all three masses of the triangle are non-zero. Consider a physical triple cut in an 
amplitude where all three corners are non-null 




Figure 2: The triple cut of an amplitude 

As the momentum invariants, A'l = /cm + ^m+i + ■ ■ ■ + etc, are all non-null, 
there exist kinematic regimes is which the integration has non- vanishing support for 
real loop momentum. 

In this section we present the contributions of various canonical forms to the co- 
efficient of the three-mass triangle integral function. As before we build our canonical 
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forms from a simple starting point. Consider, 



,,,„.iM.J^^. ,5,, 

We chose to manipulate this as if it were a covariant integral. This means effectively 
replacing, 

J dUPS J (5.4) 

then evaluating the covariant integral only keeping the coefficient of the three-mass 
triangle in the result. 

The covariant integral is a linear box function. Evaluating this we find the 
following contribution to the three- mass triangle coefficient: 



(5.5) 

2{a\KiK2\a) 

As in eqn. ( p.23|) , we could write the denominator as, 

{a\K^K2\a) = {a K^)[K^ K2]{K2a) , (5.6) 

where Ki are the null linear combinations of Ki and K2. These linear combinations 
involve irrational coefficients as in eqn. (|3.24|) . The case where (aKi) = must be 
treated as a special case. This does not arise when a denotes an external momentum 
but may arise in the derivations of more complicated canonical forms. When a = K^, 
we have. 

More complicated forms can readily be generated from this simple starting point. 
A summary of these canonical forms is given in appendix 0. As an example of the 
use of these canonical forms, the general expression for a NMHV three-mass triangle 
coefficient is given in the next section. 



6. Example: n-point Three-Mass Triangles for J\f = 1 and 
Scalar Contributions to NMHV Amplitudes 

As an example of using the canonical forms for the three mass triangle let us evaluate 
the general form for the three-mass triangle for NMHV amplitudes. The J\f = 1 
contribution was previously presented in [9]. 

For both the chiral A/" = 1 and scalar contributions, the only three-mass triangles 
which appear in the NMHV amplitude have exactly one negative helicity on each 
corner. Consider such an integral function with the following labelling: 
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Figure 3: The three mass triangle integral function appearing in the NMHV amplitudes 
indicating the labeling of the legs. 



First we take the A/" = 1 case. Evaluating this using the triple cut of fig. |] the 
triple cut integrand is, 

A(— £q, ■ ■ ■ , nil ) ■ ■ ■ ) ^1 ^ ■ ■ ■ ) ^2 i' ' ' ^^2^) ^ ^(~^2? ■ ■ ■ 5 ""^3 ? ■ ■ ■ ) ^0 '^)' 

h 

(6.1) 

where = (-q — Ki etc. The summation is over the M = 1 chiral multiplet. The 
effect of this summation is to give the scalar contribution times a p-factor, 

A(-£0,---,mr,---,£?)xA(-£?,---,m2-,---,£°)xA(-£0,---,m3-,---,£°)xp, (6.2) 

where, 



((mi 4) (m2 £2) (ma 4) - K fp) (ma fi) (mg £2))^ 
^ (mi ^l) (ma 4) ("^3 4) ("^1 4) ("^2 4) ("^3 4) 

(^oX)2 



(6.3) 



[4 4]2(mi £1) (m2 £2) (ms 4) (mi 4) (^2 4) (mg £2) ' 

and, 

\X) = \mi){m^\K^K2\m2) + Img) (mi|;riir2|m2) . (6.4) 

The cut is then, 

(mi 4) (mi £1) ("^2 4) ("^2 4) 

X 



(4 /l) (/l ■ ■ ■ «l) («1 ^1) (^1 4) (^1 /2) (/2 ■ ■ ■ W2) (W2 £2) (4 4) 

(ma £2) (ma 4) (4^)^ 



X 



{^2h){h■■■u,){u,^,){^,^2) [44]^ 

(mi£o)(mifi) (m2fi)(m2£2) {m^^2){m^^o) {X ip)^ 

{iofl){uiil) (4/2)(/24) (4/3)(n34) (414414)' 

(6.5) 

where, 

^ {Ulf2){u2h){u^h) , . 
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This can be turned into a function of io only: 

1 



(6.7) 



where, 



S = {\f,),\us),K,K2\f2),K,K2\u,),K,K2\h),K,K,\u2)}, 
Ti = {\mi),\ms),KsK2\mi),KsK2\m2),KiK2\m2),KiK2\ms),\X),\X)}, 

and we have used, 

(ha) _ {iolhhla) _ {io\K3K2\a) (£2 a) _ {io\K,K2\a) 
(iib) {io\i2ei\b) {£o\K,K2\b)' {£2b) {£o\K^K2\b)^ 



{61 



(6.9) 



This is precisely the canonical form J'^ with n = 6 as defined in appendix So the 
three mass triangle coefficient is precisely, 

bl^ = Cox4{S;T,;K,). (6.10) 

This general expression simplifies in many cases: if the rrii coincide with any of the 
Ui or fi, the Jg function simplifies to J° with n < 6. 

The scalar case is more complicated. The cut integrand is, 

(mi4)'(mi^i)2 {m2iiy{m2i2y 

X 



(4 /l) (/l ■ ■ ■ Wl) {Ul il) (£1 4) (4 /2) (/2 • • • «2) (m2 4) {i2 il) 

(ma £2)^(^3 4)' 



X 



(^2/3)(/3---W3)(%4)(4^2) 

(mi £o)'(mi (m2 £i)2(m2 £2)' {m^ £2)^(^3 4)' [^1 ^2? 



{iofl){Ulil) (^l/2)(/2 4) (^2/3)(^^3 4) (41^1^214) 

{m,ioy{msiof {m,i,f {m^i2f (m2 £2)^[^i 41^(^2 fi)^ 

(^0/l)(«l4) (^l/2)(n2^l) (/3^2)(^2/3) {^K2K^%) 

n|,)eT,(^ol/) (m2£2)'[^i^2]'(m2£i)2 



n|.)e5(^oa:) (4|i^2if3|4) ' 

(6.11) 

where, 

T2 = {|mi),|mi),|m3),|m3),i^iir2|m3), 7^17^21^3), ir3^2|mi),ir3^2|mi)}. (6.12) 
Splitting the pole and defining T' = T2 — {|mi), |mi), |r?7,3)}, the cut integrand be- 
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comes, 




(6.13) 



+ 2 (m2 1 -fCi -f^s 1 777,2) Ji (x; 7712 , mi , mi , ma ; 1 ?7T.2 ) ) 
+ J^(x; 1712, 1712, mi, mi, ma; K2\m2), K2\m2)) ) , 



where, 



Cx — 



Il\y)eT'{y^) 



(6.14) 



H\z)es-{\x)}i^ ^) 



This is a general, exphcit, rational form of the 7i-point NMHV three-mass triangle 
coefficient. 

7. Summary 

We have presented an implementation of the Unitarity method and applied it to 
the computation of the seven-point one-loop Af — 1 amplitudes. Once the canon- 
ical forms are derived, the method is algebraic. In many ways our canonical form 
approach is equivalent to alternate methods however it naturally applies to trees 
written in the spinor helicity formalism and produces results which are manifestly 
rational. 

Acknowledgements: This research was supported by the Science and Technology 
Facilities Council of the UK. 

A. Integral Functions 

The scalar box integral is. 




P 



|2 ip - Kif ip-Ki- K2f {p + X4) 



1 



2 ' 



(A.l) 



where Ki is the sum of the momenta of the legs attached to the i-th corner. If a 
single leg is attached then is null. The form of the integral depends upon the 



number of the Ki which are non-null, ^ 0. We often misname these massive legs. 
The integrals are functions of the non-zero Kf and the invariants, 

S={K^ + K2)\ T = {K2 + K,y. (A.2) 

For convenience, we always define these integrals taking leg 1 as massless and leg 4 
as massive. 

It is convenient to define the scalar box function, F, 
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and the symmetric 4x4 matrix S has components {i, j are mod4). 
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In terms of these variables, the relationships between the scalar box functions and 
scalar box integrals are given by, 
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With the labelling of legs shown in fig. the scalar box functions expanded to 
(9(e°) in the different cases reduce to, 
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Figure 4: The different types of scalar box. 
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Our seven-point expressions do not need the explicit form of the four mass scalar 
box. 

At this point we must discuss a suitable basis for expressing the amplitudes. We 
could use the basis ( p. . 1|) , however this is not the most efficient option. By choosing 
a suitable basis of box functions we can considerably simplify the structure of the 
triangle coefficients. 

Triangle integral functions may have one, two or three massless legs: /g™, /|™, 
/|™. The one-mass triangle function depends only on the momentum invariant of 
the massive leg Ki and is, 



(A.ll) 



while the two-mass triangle function with non-null momenta Ki and K2 is. 
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Both of these integral functions contain ln(fC^)/e IR singularities. The key point is 
that the IR singularities of the amplitudes must be [39] , 

j^N=l chiral _ ^^tree ^[0] _ ^^tree j^oN 

so the ln(A'^)/e singularities must cancel. This constraint effectively determines the 
coefficients of and /|™ in terms of the box coefficients. 

Specifically the one- and two-mass triangles are linear combinations of the set of 
functions, 



G{-K')^rr ' / , (A.14) 



with, 



/r=G(-iffl , 

The G{—K'^) are labeled by the independent momentum invariants and in fact 
form an independent basis of functions, unlike the one and two-mass triangles which 
are not all independent. 

In practice we need never calculate the coefficients of the G functions once we 
know the box coefficients. The only functions containing ln(s)/e terms are the box 
functions and I^"^ and /|™ so, 

J2ajX^K2y, + bG^^^^^0. (A.16) 

This equation fixes the single be in terms of the a^. The simplest approach to 
implement this simplification is to express the amplitudes in terms of truncated, 
finite ^-functions [3,40,23]: 
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The J\f — 1 and scalar amplitudes can then be expressed as, 

^i-ioop ^ J2 a^r + ^ \?;^fp' + ^ Ckll + R, (A. 18) 

with no /g"* or 7|"* present. 
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In transferring to this set of basis functions, the coefficients of the JF*, ojf, are 
simply related to the coefficients of the scalar box integrals, a/, by, 
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The remaining integral functions are the three-mass triangles /I™" as given, for 
example, in ref. [41] and the two-point bubble function, 

UP) = 7 + 2 - In(-P2) . (A.20) 
B. Canonical Forms for Triple Cuts 

In this appendix the contributions of various canonical forms to the three-mass tri- 
angle coefficient are given. The triple cut is labelled as in figure p|. 
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r,^ii _ (6|[iri,ir2]|rf)(c|[iri,ir2]|a) + A3(&rf)(ca) (c;6)(rfc)(a|[iri, iT^ 

J, [d, a, D, c, lA.lJ - 2A3(rf|iriir2M) 2miKM? 

(B.l) 

where. 

As = ^{KrK^f-^KlKl = {Klf+{Klf+{Klf-2{KlKl+KlKl+KlKl). (B.2) 

The above expressions are valid for £ = ^o? ^ = and i = £2 since, 
{a\[K,,K,]\b) {a\[K^,K2]\b) {a\[K2, Ks]\b) 



2{a\KiK3\a) 2{a\KiK2\a) 2{a\K2Ks\a) 
We usual, extend these to, for example. 



(B.3) 



The J° can be expressed in terms of J° as, 

n 

J'^[{d,}- {a, b, q}; {A^,}] = C^J'Ad^■, a, b, c„; {K,}], (B.5) 

i=l 
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with, 

lYiZlicjdi) 



We also have terms which are of order 

£^ = [A\£o\a) E'o[a;A;{Kj}] = -[A\ao\a) , 
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' {tod) 
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(B.7) 

where, 

+ Kl - i^l) + Kl - ) 
ao = ^- + ^- i^a • (B.8) 

Expressions for £ — ii are obtained by replacing oq by oi where, 

= —K^ + \ —Ki^ = Oo - ■ (B.9) 

The i = £2 expressions are obtained in a similar fashion. 
Finally for order we have, 

[A\eo\b){ioc){eod){eoe) 

{eo\K,K2\£o){iof) ~" 
TMfh ^ . ^ [A\ao\b){f\[K„K2]\c){d\[K^,K2]\e) 
J, [f ; b, c, d, e, A, |i^,}J - 2^ 4As{f\K^K2\f) 

^ [A\ao\f){b\[K,,K2]\c){d\[K,,K2]\e) 
^ 2AAs{f\K,K2\f) 

[A\ao\f){fb){fc){d\[K,,K2]\e) 



^ 12{f\K,K2\f)^ 
^ [A\ao\b){f\[K,,K2]\c){f\[K,,K2]\d){f\[K,,K2]\e) 
^ AAMK1K2W 
^[A\ao\f){fb){fc){f\[K,,K2]\d){f\[K^,K2]\e) 



m\K,K2\fr 



(B.IO) 



- 29 - 



For QCD amplitude we generically need forms of order 
£^^[A\io\a)[B\io\b)^ 
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Cubic and higher order terms can also be evaluated. For example, 

[A|£o|a)[B|£o|&)[CKo|c) ^-i5^[^|ao|a)[B|ao|6)[C|ao|c) 



(B.12) 

however we can, in general, avoid these in QCD calculations. 



C. AU-n Expressions for Coefficients 
C.l Box Coefficients 

Many of the box coefficients that appear in the seven-point M — \ amplitudes are 
special cases of general n-point expressions. We gather these together here and give 
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their specialisations to the box coefficients of section ^. We denote the external legs 
contributing to the non-null momentum Ki by {/j, ■ ■ ■ , Mj}. 

The three-mass boxes in a NMHV M = 1 amplitude vanish unless exactly one 
negative helicity gluon is attached to each non-null vertex. We denote this single 
negative helicity gluon on vertex i by 
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where, 
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'H = {nil 1112) {d\K2K3\m3) + ("^3 "^2)('^|-f^4-f^3|"^i)- 



(C.2) 



These box coefficients are close in form to those of AT = 4 Yang- Mills [13, 15]. 

All of the NMHV two mass boxes we need can be found using Generalised Uni- 
tarity taking MHV and MHV tree amplitudes as the inputs. 
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With the labelhngs given in the figure, the coefficients are, 
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(C.3) 

where, P = + h = + /cmj. 

One mass boxes in A/" = 1 amphtudes have either one or two negative hehcity 
external legs attached to the three-point corners. In the latter case, the third negative 
helicity leg of any NMHV amplitude attaches to an MHV corner and a general form 
exists: 





^Im [b\Ki\m3)'^[mi\Ki\m3)[m2\Ki\m3)Srn,bSbm2 

° 2Kj[mi m2f[mi\K4\m2 + l)[m2\Ki\ui) YTi^jli^^ + 1) 



When there are two negative helicity external legs attached to the massive corner, 
NMHV tree amplitudes are required and we must consider each helicity configuration 
separately: 
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(C.5) 

C.2 Particular Bubble Coefficients 

Our first example is the bubble coefficient in an jV" = 1 MHV amplitude. Coefficients 
of ln(— P^...j) vanish unless exactly one of the two negative hehcity legs hes within 
Pa-.-b- The non-vanishing coefficient of ln(— P^...j,) is then, 

(mi m2)^(&6 + l)(a — 1 a) n n \ 
X Hi[a - 1, a, 6, 6 + 1; mi, mi, m2, mg; Pa-bJ • (C.6) 

The second example is relevant for the seven-point NMHV amphtudes. Consider 
the case where one side of the cut is MHV and the other is MHV. Let m be the 
single negative helicity on the MHV side which contains legs a - ■ - b and p be the 
single positive helicity on the MHV side. The coefficient of ln(— P^...^) is then just, 

[p\Pa...b\mf 

H4[a,b,Pa...b\a - 1], P„...b|& -M]; m, m, P„...b|p], P„...6|p]; P„...6] . 
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Since many of the di coefficients in the seven-point NMHV amphtudes are special 
cases of this generic expression specialised to seven-point we define, 

Co[a,b,c,d,eJ,g;p,m] = 
\v\P Itti)'^ 

{cd){de){e7)[9a][ab]P^^, ^ ^' ^-^1^1' ^-^1^1; ^' ^' PgaM Pa-.] ■ 

(C.8) 

D. Seven-point Tree Expressions 

We have: 

A : A(s,, S2, 3-, 4-, 5+, 6+, 7+) = T,i + T^, + + , (D.l) 
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B : A{s,, S2, 3-, 4+, 5", 6+, 7+) = T,^ + Tf, + T^^ + + 7f + Tf , (D.3) 
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C : A{s,, S2, 3-, 4+, 5+, 6", 7+) = Tg + Tg + Tg + + T§ + Tf , (D.5) 

with, 
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D : S2, 3-, 4+, 5+, 6+, 7-) = T/^ + + T^^ + T^^ + , (D.7) 

with, 

[6|Pn|3)2(32)[16P ^ [6|P7l|3)^^'^ 
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E : S2, 3+, 4-, 5-, 6+, 7+) = T^^ + T,^ + Tf^ + + 7f + Tf , (D.9) 
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with, 
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P : A7(si, S2, 3+, 4-, 5+, 6", 7+) = Tf, + T[, + T^^ + + Tg^ + Tf , (D.ll) 

with, 
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' (12)(23)(45)(56)[7|Pi23|3)[7|Pi23|4)(l|P23P45|6)tl23t456 V [71^2312)/ ' 
^ (1 3)2+2^(2 3)2-2ft ^ ■ 

(D.12) 



E. Functions 



The function Da was given in section [4.5| . The remaining five Dx functions are given 
by: 
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DB[a,b,c,d,e,f,g] = 
1 



G^ia, g, PabPcd\e) , Pab\c], Pefg\d];b, e, e, X^ia, X^ia, Pab\d];d; Pefg;Pab] 

^(x yu U/j 

+ J^^^f^^He[(^,9,Pf9\e],Pefg\d],Pf^^^^ 

H5[a, g, PabPcde\f),PabPcde\e), PfgPcde\c);b, a, /, PfgPcde\b),PfgPcde\b);Pab] 
H5[a, g, PabPcdelf), PabPcde\e), PfgPcde\c)\h, b, /, PfgPcde\b),PfgPcde\b);Pab] 
H^ia, g, PabPcdelf), PabPcde\e), PfgPcde\c);b, a,g, PfgPcde\b), PfgPcde\b);Pab] 
H5[a, g, PabPcdelf), PabPcde\e),PfgPcde\c);b, b, g, PfgPcde\b),PfgPcde\b);Pab] 



{ef){fg){ab)[cd] 
\d 

{ce)naf] 
{ab){cd){de){f g)tcde 

{ab){cd){de){f g)tcde 
{ce)^[ag] 

{ab){cd){de){f g)tcde 
{ce)^[bg] 



{ab){cd){de){f g)tcde 

{ce)^[ga]^{a b)^ 

H2 [a, PabcPdelf); b, c; Pab\ 



{a b) {d e) (e /) [g\Pabc\c) [g\Pabc\d)tabctdef 

[f9]Hbc)^ 



{ab){cd)[ef]tefg[g\Pef\d) 
where, 



H2[a,Pfg\e\;b, c;Pab] , 

(E.l) 

\XBia) = [da]{be)\a) + [d\Pbcd\e)\b) . (E.2) 



Dc[a,b,c, d,e,f, g] 
1 



[bc]{de){ef){ga 

[cdf{ba)^ 



GbiPgalb], f, Pdef\c], PgaPbc\d) , Q', 6, 6, Xcia, ^cia, O, -Pga|c]; C; Pbc, P, 



ga\ 



{ef){ga) 
1 



Half, g, PgabPcd\e),PefPcd\b), Pdef\d], Pdef\c]; a, b, e, e, e, e; P, 



{ef){bc){cd){ga)tbcd 

'XG5[f, g, PefPcd\b) , PgaPbcd]^) , PgaPbcd\d); a, 6, 6, Xcic, Xqic, PgaPbcd\b); Pbcd\b\, Pef', Pgo\ 

{he)\ag)^[gff ^ r p I \ f Pi 
H2[a, Pbcd\e);f, g; Pga] 



{b c) (c d){de) {g a) [f \ P^de \ b) t fgaSga 
^ [f\Pcd\e)'{ba)' 

+ / ,\/, \r.|D I \r.|D |,\, : 7 rH2[g, PgabPcd\e); O, b] Pga\ 



-H,[g,P,f\d\,aMPga], 



{bc)[de][ef][f\Pde\c){ga)tdef 

(E.3) 
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where, 



\Xcia) = -{ea){[cg]\g) + [ca]\a)) + [ch]{he)\a) , 
l^cic) = (e6)t;,cd|a) + (ea)((6|Pcd|^]|5) + {h\Pj\a]\a)) . 



Dnia, b,c,d,e,f,g] = 

{ab){cd){de) [f g] ^^blg], PabPfg\e), Pcde\f]', Ct, C, Xdi, Xdi, Pab\f]', /; -Pcde; Pab] 

■I P , 1^3 

-H2[b,Pcde\f];a,g;Pab] 



{ag)V\Pcde\ 


|c)3 


{a b){cd){de){e\ 


PcdPab \ g)tcdetgab 



(a6)(e/)(/(?)[cd][c|P„fe|^)(^|P„bPed|e 
1 

+ 



{ab) {ef){fg){cd)tefg 

X G^lb, PabPefgl^^), PabPefg\d) , PcdPefld)', (^1 C, PabPe f\a),Xu^,Xoi\ Peflg) ) Pcdi Pab\ 

(E.5) 

where, 



\Xni)^-\a)[f\Pga\c) + \b)[bf]{ac), 
\Xn4) = |a)[d|Pe/b)(cd) + \c)[b\Pef\g){ab) . 
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DE[a,b,c,d,e,f,g] - 
1 

{cd){ef){fg)sab 



X 



/ [f\Peb\dy 


G5 


a,g, 


C^Pfgk 


,Pcd\e] 


; d, 


e, b, /, 6, / 


a; Pefg 


Pab]\ 


+[g\Peb\dr 




a, 9 


c,Pfg\e 


,Pcd\e\ 


■d, 


e,b,g,b,g; 


a; Pefg-, 


Pab] 


+{[fa]{bd)r 




a, 9: 




,Pcd\e] 


d, 


e,a,f,a,f 


; a; Pefg 


■,Pab] 


+ {[ga]{bd)y 


G, 


a, 9, 


c,Pfg\e 


,Pcd\e\ 


d, 


e,a,g,a,g 


a-', Pefg 


Pab] 


+2[m,\d)[g\P,,\d) 


G, 


a,g 


c,Pfg\e 


,Pcd\e] 


; d, 


e,&, /, b,g; 


'^'-1 Pefg'-i 


Pab] 


+2[m,\d)[fa]{bd) 


G, 


a,g. 


c.PfM 


,Pcd\e] 


d, 


e, b, /, a, / 


a.; Pefg 


Pab] 


+2[m,\d)[ga]{bd) 


G, 


a,g, 


c,Pfg\e 


,Pcd\e] 


; d, 


e,b, f,a,g 


0-', Pefg 


Pab] 


+2[g\Peb\d)[fa]{bd) 


G, 


a,g, 


c,Pfg\e 


,Pcd\e] 


d, 


e,b,g,a,f 


a; Pefg 


Pab] 


+2[g\Peb\d)[ga]{bd) 


G, 


a, 9, 




,Pcd\e] 


; d, 


e,b,g,a,g 


a; Pefg 


Pab] 


\+2[fa]{bd)[ga]{bd) 


G, 


a, 9, 


c,Pfg\e, 


:Pcd\e\ 


d, 


e,a,f,a,g 


a; Pefg 


Pab] J 



{ed) 



{cd){ef){fg){df)sab 

( {mbW 

^{\9\Peb\d)f 

+ ([/a](6c^))2 
^{\9o\kbd)f 



X 



G^\a, g, c, Pfg\e],Pcd\e]] d, f, b, f, b, f 
Ggfa, g, c, Pfg\e\, Pcd\e]; d, f, b, g, b, g 
G5[a, g, c, Pfg\e],Pcd\e];d, a, /, /, a, / 
G^la, g, c, Pfg\e], Pcd\e]; d, f, a, g, a, g 
+2[f\P,,\d) [g\Peb\d) G,[a, g, c, Pfg\e], PM;d, f, b, f, h, g 
+2[f\P,b\d)[f a] {bd) G,[a, g, c, Pfg\e], Pje]; d, /, b, f, a, f 
+2[f\Peb\d)[g a] (bd) Ggfa, g, c, Pfg\e], Ped\e]; d, f, b, f, a, g 
+2[g\Peb\d) [f a] {bd) G^la, g, c, Pfg\e], Pcd\e]; d, /, 6, g, a, / 
+2[g\Peb\d)[g a\{b d) G^[a, g, c, P/g|e], Pcd\e\;d, /, 6, g, a, g 
\+2[fa]{b d) [g a] {b d) G^ia, g,c, Pfg\e], Pcd]e];dJ,aJ,a, g 



PfgXEs 


Pab]\ 


PfgXEz^ 


Pab] 


PfgXE3 


■,Pab] 


PfgXEz 


Pab] 


PfgXE:v 


Pab] 


PfgXEz 


Pab] 


PfgXEi 


Pab] 


PfgXEs 


Pab] 


PfgXEs 


Pab] 


PfgXE3 


Pab]) 



+ 



+ 



{ab){ef){fg){df)[c\Pae\f) 





PedeW 


{fg){ab)[cd][de][c 


Pde\f)tcde 



{de^lga] 



(c d) {d e) (e /) [a b] [g\Pab\c)tgab 

{def[g\Pdef\by 
{a b){ef)[g\ Pdef \c)[g\Pef\d) tabctde / 
[f9?{bdf 



G,[g 1 ^1 PabPcde I/) ) PabcPde 

\f)]b, f, f, Xeu, ^Elb] C; PfgXEz'i Pab\ 

H2.[a, g, Pcd\e];b, Pde\c], Pde\c]; Pab] 

H2[a,PabPgab\f);b,Pab\g];Pab] 

Hala, C, PabcPde\f);b, Pdef\g], Pdef\g]; Pab] 



{ab){cd)[ef][g\Pef\d)t 



efg 



H3[a,c,Pfg\e];b,d,d;Pab] , 



(E.7) 
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where, 

\XEit) = [cg]{gf)\b) + [ca]{ab)\f), 



l^x.3-^l/)(e|. 



, ,u ,}/ g, Pefg\d\,PabPcd\e);f, f, a, Xpia, Xpia, q Pcd; Pab] 

{ab)[cd]{ef){fg) 
[ec] 



{fg)[de]{ab) 



H7[b, g, C, Pfg\e\, PfgPde\c),Pabc\d], PabcPde\f);a, C, /, /, PabM, Xpu,, Xfu,; Pab] 



+ ^yg^[Je](a6) "^^ ^' "^-^^ ^'^'^ ' ' '^"'"^'^'^^ ' ' ^"^^1^1 ' ^f9 1^] ; /' /' Pabc\e] , Pab\e] , Xpu, Xpu; Pab] 
1 



if g){ab){cd){de)tcde 

X G5[b,g,PfgPde\c) ) PabPcd 1 6) ) PabPcde 

\fy,a,f,f 

) PabPcde \d) , Xpic, Xpic, Pcde\d]', Pfg Pab] 

{df)'[gb]'{ba) 



(c d) {de) (e /) b|P(ie/ \c)Sabtgab 

{df)Hg\Pabc\a]^ 



{ab){de){ef)[g\ Pahc \c)[g\ Pabc \ d) tabcUef 

[ge]Hda)^ 



H2[Pab\b),Pcde\f);g,Pab\a);Pab] 

Hsib, C, PabcPde\f);Pabc\g],Pabc\g],a; Pab] 



I u\i fMf ir id^CaT — H2.[b,c,Pefg\e]]d,d,a\Pab], 
{a b) (c d) [e /] [/ g] [g\Pefg\d)tefg 

(E.9) 

where, 

|^Fia) = |/>[c6](6a) + |a)[c|Pe/p|/), 

\XF,b) = \f)[ed]{da) + (/a)(|/)[e/] + \g)[eg]) , (E.IO) 
I^Fic) = \a){fg)[g\Pcde\d) + |/)(a6)[6|Pc<ieM) ■ 
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